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Hybridization and multipole orders of 4/ electrons in Pr skutterudites 

Yoshio KuRAMOTO*, Annamaria Kiss, Junya Otsuki and Hiroaki Kusunose 
Department of Physics, Tohoku University, Sendai 980-8578 

Characteristics of hybridization and multipole orders of 4/-electrons in Pr skutterudites are explained 
in terms of a pseudo-quartet composed of crystalline electric field (CEF) singlet and the triplet. It is 
shown that the contrasting behaviors observed in PrFe4Pi2 and PrOs4Sbi2 are ascribed to the difference in 
triplet wave functions and the CEF splittings. Since a macroscopic degeneracy remains even in the ordered 
phase with the Fj-type antiferro-quadrupole (AFQ) order, a model with strong AFQ fluctuation and static 
monopole/hexadecapole order is proposed for PrFe4Pi2. The remaining cubic symmetry explains qualita- 
tively the behavior of staggered magnetization observed by NMR and neutron scattering. For identification 
of possible hexadecapole orders in PrFe4P|2 and PrRu4Pi2, eightfold intensity pattern is predicted in the 
azimuthal angle scan of resonant X-ray scattering. In PrOs4Sbi2, the ferromagnetic exchange coupling 
with the conduction band does not lead to magnetic Kondo effect, but a momentum-dependent quadrupole 
coupling can give rise to enhanced effective mass, which should be sensitive to disorder. 
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1. Introduction 

There are puzzling properties in filled Pr skutterudite com- 
pounds. In this paper we concern with mechanism of the fol- 
lowing behaviors: PrFe4Pi2 undergoes a phase transition at 7 
K in zero field. ' The susceptibility shows a cusp at the tran- 
sition, but there is no Neel order in the ordered phase. Then 
an antiferro-quadrupole (AFQ) order has been considered as a 
plausible explanation.^ Indeed a staggered magnetic moment 
and the lattice distortion have been found by neutron diffrac- 
tion.^ However, in contrast with the cas of PrOs4Sbi2,"^ there 
is no component perpendicular to the field of the staggered 
moment.^- ^ If the order parameter is AFQ, one may naturally 
expect finite perpendicular component in some directions of 
the field. Thus proper identification of the order parameter 
remains to be done. On the other hand, PrOs4Sbi2 has well- 
defined CEF excitations as seen by neutron scattering,'*'^ and 
by the Schottky-type anomaly in the specific heat.^'** Near the 
Schottky peak around 2 K, superconductivity sets in with a 
huge discontinuity in the specific heat. The problem here is 
the origin of the Cooper pair. More modestly, one may ask 
the mechanism which produces heavy electrons forming the 
Cooper pair. These two Pr skutterudites also show contrast- 
ing behaviors in the temperature dependence of the resistivity; 
PrFe4Pi2 shows a typical Kondo effect, while PrOs4Sbi2 only 
a mild shoulder in the temperature dependence.^ The inelas- 
tic neutron spectrum of PrFe4Pi2 is characterized by broad 
quasi-elastic features like other heavy-fermion systems, with- 
out well defined CEF excitations. 

In this paper we propose a possible route toward under- 
standing these fundamental issues. We begin with brief review 
of our previous results which explain the contrasting behav- 
iors in terms of the character of CEF states. Then we pro- 
ceed to rather speculative new ideas. These include proposal 
of a novel multipole order in PrFe4Pi2 in which fluctuating 
quadrupoles are not frozen even in the ordered phase. The 
frozen multipole is a hexadecapole which should be probed by 
resonant X-ray scattering. In addition, for PrOs4Sbi2, we pro- 
pose that the induced quadrupolar fluctuation with momen- 
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tum dependent coupling with conduction band is responsible 
for the heavy mass. 

2. Local electronic states and hybridization effects 

2.1 CEF splitting by p-f hybridization 

Filled skutterudites RT4X12 form a bcc structure with the 
space group Im3. The R site has the local symmetry Tf, which 
has no four-fold rotation axis." In this symmetry, the 4/^ 
Hund's-rule ground states ^7/4 of Pr"*^ split into a singlet Fi, a 
non-magnetic doublet r23, and two r4 triplets. Of these, two 
r4 triplets are written as F^' ' and F^^', which are linear combi- 
nations of triplets r4 and Fs in Oh- Under this crystal symme- 
try, the CEF potential is written as 
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in the standard notation. ' ' The term yO'^ mixes the F4 and F5 
triplet states in the point group Oh- 

There are two main sources for the CEF splitting: the 
Coulomb potential from ligands, and hybridization between 
4/ electrons and ligands. In the point charge model, CEF co- 
efficients A4, and A^ in eq.(l) are determined by coordina- 
tion of the charge and the radial extension of the 4/ wave 
function. Explicit results have been obtained in ref. 14 tak- 
ing the effective charge Zt of a transition ion as the param- 
eter and requiring the charge neutrality with trivalent Pr To 
estimate the CEF potential, we use the lattice parameter in 
PrOs4Sbi2,'^ which gives the Pr-T distance dx — 4.03A with 
T=Sb, the Pr-X distance c/p = 3. 48 A with X=Os and the X- 
Pr-X vertical angle 29q - 49.2°. These data will be used to 
derive the results presented in Fig. 1 . 

Another important mechanism for CEF splittings is the co- 
valency effect, or hybridization between localized and lig- 
and orbitals. According to band calculation,'^ the conduction 
band striding the Ferrrd level is formed mostly by the molecu- 
lar orbital a,, formed by 12 pnictogens around each Pr. The 
other relevant molecular orbital f„ form two energy bands 
around a few eV above the Fermi level, and three a few eV 
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below. We neglect contributions from tu bands to CEF split- 
tings because of the larger excitation energy. Therefore we 
first consider hybridization of the form 

H^yb = VaYjflc^+^-C-' (2) 

where c,j- annihilates a conduction electron in the Wannier or- 
bital with the a„ symmetry at the origin, and fl creates an 
4/ electron with the same (a„) orbital symmetry. We take the 
hybridization parameter Va real. Hereafter we adopt the Mul- 
hken notation such as Uu for orbital symmetry, and the Bethe 
notation such as Fi for a double-group representation with 
spin-orbit coupling. In the second-order perturbation theory, 
the eff"ective interaction is given by 

1 

Hefi = PHbyh ^ QHhybP, (3) 

where P is the projection operator onto 4/^ states, and Q = 
1 - P. We first deal with such part of the second-order hy- 
bridization that is diagonal with respect to the conduction 
states. Diagonalization of H^g with this constraint gives the 
CEF wave functions and their energies. We assume the fol- 
lowing intermediate states: (i) 4/' and an extra electron in 
vacant states, and (ii) 4/^ and extra hole in filled states. For 
simplicity, we neglect the multiplet splittings in both kinds 
of intermediate states. We note that 4/' and 4/^ intermediate 
states give opposite contributions to the level splitting. Hence 
the sequence of CEF levels is determined by competition be- 
tween both intermediate states. The hybridization is parame- 
terized by the Slater- Koster parameters (p/cr) and (pfn). It 
turns out hybridization with a„ comes only from (pfn). 

We combine both contributions to CEF splittings by point- 
charge interaction and by hybridization. To a good approxi- 
mation we can assume the half-fiUed a„ band. We introduce a 
parameter 

1/A_ = I/A3 - 1/Ai. 

Note that 1/A_ = in the case of A3 = Ai, and 1/A_ be- 
comes negative if 4/^ is the dominant intermediate state. For 
the point charge parameters we tentatively take Z, = 2. Fig- 
ure 1 shows computed results as a function of strength of hy- 
bridization, (pf7i)^/A-. It should be emphasized that the the 
closely lying singlet and triplet are realized only with simul- 
taneous action of point-charge and hybridization interactions. 

The lowest triplet Ft, which is either F^^' or F^'^', can be 
written in terms of Oh triplets F4 and F5 as 

|Ft,m) = Vw|F4,ot) -I- Vl -w|r5,OT), (4) 

where m = ±,0 specifies a component, and w with < 
w < 1 gives the weight of r4 states. We have chosen the 
phase of wave functions so that positive coefiicients give the 
lower triplet. The level repulsion between F^ ^ and F^^^ around 
(pfnf/A- = llOKis due to mixing of wave functions, and 
is a characteristic feature in the point group T/,. Possible lo- 
cation of effective hybridization in some representative Pr 
skutterudites is schematically shown by the lines (a)-(e) in 
Fig.l. In the case of PrFe4Pi2, the lattice constant is 7.81A as 
compared with 9.3A of PrOs4Sbi2. Hence (pfn) in PrFe4Pi2 
should be larger. At first sight, the larger hybridization in 
PrFe4Pi2 seems incompatible with this leftward shift (a) from 
PrOs4Sbi2 around (d). However, since the smaller lattice con- 
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Fig. 1. CEF level structures derived from hybridization and point charge 
potential as a function of (/)/;r)^/A_. The level sequence qualitatively cor- 
responds to: (a) PrFe4Pi2; (b) PrRu4Pi2 in the high-temperature phase; (c) 
Prl site in PrRu4Pi2 in the low-temperature phase; (d) PrOs4Sbi2; (e) Pr2 
site in PrRu4Pi2 in the low-temperature phase. See text for details. 

Slant in PrFe4Pi2 favors the 4 /' excited states rather than the 
4/^ states, the parameter Ai can become smaller, leading to 
negative A_. Since the hybridization with a„ orbitals alone 
cannot take into account the parameter y in eq.(l), the si- 
multaneous crossing of three levels near (a) should be taken 
with caution. In other words, only the Coulomb potential is 
the source of the y parameter in our result. From the study of 
phase diagram in magnetic field,'^ it seems that the level F23 
is situated higher in this region. As a result, if high pressure 
drives the system toward further left in the Fig.l, it is likely 
that the triplet becomes the ground state instead of the dou- 
blet. 

It is reasonable that the larger cage in PrOs4Sbi2 prefers 
4/^ intermediate states. In PrOs4Sbi2, the triplet forming 
pseudo-quartet is almost of the Fs-type as shown by (d). On 
the other hand, the CEF levels in the high-temperature phase 
of PrRu4Pi2 corresponds to (b), which bifurcates into (c) and 
(e) in the low-temperature phase. Namely in one sublattice 
called Prl, the triplet changes the character shghtly, while in 
the other sublattice called Pr2, the triplet becomes the ground 
state. Below the phase transition, the system behaves like an 
insulator. It may be possible to describe the transition as the 
antiferro-hexadecapole order as discussed elsewhere. '^'^ 

2.2 Exchange interactions with conduction electrons 

We are interested in the matrix element between the singlet 
and the triplet in the T/, group. As we have seen the triplet is 
a linear combination of F4 and F5 in the Oh group. Hence it 
is convenient to analyze the symmetry properties of triplets in 
the cubic case. The direct products of pseudo-quartets can be 
decomposed as 

(Fi e r4) ® (Fi e r4) = Fi -h 2F4 -n (f, + r3 + F4 -h Fs), (5) 

(Fi e Fs) ® (F, e Fs) = Fi + 2F5 + (F, + Fj + F4 + Fj), (6) 

where F4 has the same symmetry as the magnetic moment. 
The product of Fi ® F4 produces newly two F4 representations, 
of which one is time-reversal odd and the other is even. The 
odd representation represents the magnetic moment, while the 
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even one the hexadecapole moment. On the other hand, the 
pseudo-quartet Fi © Fs does not produce a new F4 represen- 
tation. Physically, this means that r4 as the first excited states 
gives rise to a van Vleck term in the magnetic susceptibiUty, 
while F5 does not. 

It is convenient to introduce the effective moment opera- 
tors in the pseudo-quartet in Th- In the case of w = (pure 
F5 triplet), the magnetic moment within Ft is the only relevant 
quantity. The vector operator in this case is written as Z*. On 
the other hand, another vector operator is necessary to de- 
scribe the magnetic moment of van Vleck type. The operator 
AT' act on the triplet states, and connect the singlet and the 
triplet. The effective interaction within the pseudo-quartet is 
given by 

//exc = (/t^ + kX') ■ \ '^I'^-P'^P- (7) 

where the conduction electron spin Sc is explicitly written in 
terms of fermion operators. The coupling constants are pro- 
portional to and negative. The origin of the ferromag- 
netic exchange is traced to the Hund rule involving the spin- 
orbit interaction; the dominant orbital moment is pointing op- 
positely from the spin moment. In other words, the spin ex- 
change is antiferromagnetic as is usual for a hybridization in- 
duced exchange. 

Alternatively, one can use the pseudo-spin representation 
as in ref.l9: 

where Si and S2 are spin 1/2 operators. The pseudo-spin rep- 
resentation naturally leads to commutation rules among X^ 
and Z^ The exchange interaction in terms of pseudo-spins is 
given by 

//exc =(/|Sl +/2S2)-*c, (8) 

where Ji = h + h and J 2 - h - h- It turns out that J2 be- 
comes positive for as w deviates from 0. The emergence of 
antiferromagnetic exchange is due to the particular CEP level 
structure in Pr skutterudites. J2 is almost negligible in the pure 
F5 case (w = 0), and becomes an order of magnitude larger as 
w increases toward unity, i.e., toward pure F4. It is likely that 
the Kondo-type behavior seen in PrFe4Pi2 originates from a 
F4-type triplet together with a small singlet-triplet sphtting. 
Nemerical results have been reported for dynamical proper- 
ties of the pseudo-quartet^" with use of the NCA.^' 

2.3 Quadrupole interactions with conduction electrons 

Although the conduction band is of a„ type, it can acquire 
the coupling with quadrupole moments through mixing with 
f„-type calcogen bands. We start with the most general form 
of p-f hybridization as given by 

^hyb = X X [^r(QiS)Pr(„)y^/ros)v(r + H.c] , (9) 

r(afi) vcr 

where the indices a,j8 distinguish different sets with the same 
irreducible representation F. A member in the irreducible rep- 
resentation is specified by y, and cr denotes a spin compo- 
nent. Let us consider nine Pr sites making up a body-centered 
cube each of which is surrounded by 12 pnictogens. A 4f or- 
bital state 3z^ - 5, for example, in the body centered Pr at 



R - (000) has hybridization with surrounding pnictogen or- 
bitals F = V = z, for which five different sets have been 
enumerated.^'* 

Each tu molecular orbital mixes with the eight au orbitals 
surrounding Pr comer atoms at/? = (±l,±l,±l). Sunomation 
over eight sites gives the form factor as 

F^(k) - VSicosfe^cosfe,, sinA:^, (10) 

where k is the wave number of the conduction electron. Other 
components Fy{k) with v = x,y are obtained from the case 
of V = z by cyclic rotation of k^, ky, k^. Let us introduce a set 
V - x,y,z of new Wannier orbitals each of which is created 
by the operator 




Then the effective hybridization through orbitals is given by 
^' = ^ Z Z f^^^-/^- + H.C.1 , (12) 

j3 vcr 

where yS distinguishes two different sets with the 4 symmetry 
in 4/ orbitals such as 3z^ - 5 and z{x^ - y^). The effective hy- 
bridization parameter Vp is determined by Vtqia), the intersite 
mixing between f„(a) and Uu, and their energy separation. 

Following the same line as derivation of the effective mag- 
netic exchange, the second order effects of eq.(12) gives 
quadrupole couphng. To describe the couphng we intro- 
duce the operator O^^y where (y) = (1 5) represents 
the quadrupole from the CEF singlet-triplet transition, while 
(7) = (5^5) comes from the intra-triplet transition. We 
shall argue later that the latter quadrupole operator is a key 
to understand the strange heavy-fermion behavior found in 
PrOs4Sbi2. The effective quadrupole coupling between 4/ 
and conduction electrons is given by 

y Ijy (T 

We shall evaluate these two coupling constants I^^^^^ quantita- 
tively in separate work. 

3. Fluctuating quadrupoles 

3. 1 Triplet pair with AFQ interaction 

An AFQ order of the Fa-type involving a CEF triplet in- 
volves a subtle problem concerning the degeneracy in the or- 
dered phase. The normalization of relevant operators is de- 
fined by 

Ol = 3-^'^[3J^ - J{J + 1)J, Ol^Jl- J^. (14) 

In the case of F4 triplet, nonzero matrix elements are given by 
(0|(9°|0) = 28/ V3,<+|C»^|±) = -14/ and (±10^1+) = -14. 
For identifying the degeneracy problem, which is intimately 
related to the ordering pattern in PrFe4Pi2, we take a two- site 
system where the triplet r4 is present at each site A or B. The 
interaction Hamiltonian is given by 

Hr,=gi(0l^0% + 0lj,0l,,). (15) 

We recombine three states in the triplet as 

\z) = 10), \x) = ^(1+) + I-)), \y) = ^(1+) - I-)). (16) 
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Then we take the nine-fold basis \k)A\l)B for the two sites, 
where k, I are x, y or z. By diagonalization of Hrs we find 
that the nine-fold degeneracy splits into six-fold and three- 
fold degenerate states. With antiferro-type interaction ^3 > 
0, the ground state is the six-fold multiplet with energy 
-196^3/ V3, and the three-fold degenerate excited level has 
energy 784^3/3. The six-fold degeneracy comes from the 
number of ways to choose difl'erent orbitals for |A:)a|/)b with 
k i= I. Because of the absence of ofl'-diagonal elements in the 
two-site Hamiltonian, the mean-field theory at zero tempera- 
ture becomes exact for the two-site system. The bcc lattice can 
be separated into A and B sublattices. Then with only the r3- 
type intersite interaction, the degenerate ground state of the 
bcc lattice should be given exactly by the mean-field theory. 

Thus a model with nearest-neighbor AFQ interaction of 
r3 quadrupoles has a degeneracy with respect to diff'erent 
kinds of quadrupolar patterns. In this sense, the r3 antiferro- 
quadrupolar ordering model within the triplet state is similar 
to the three-state antiferromagnetic Potts model, which pos- 
sesses a macroscopic degeneracy in the ground state. 

3.2 Interaction with conduction spins 

We now ask how the degeneracy should be broken in ac- 
tual PrFe4Pi2. Experimentally, distortion of Fe sublattice has 
been observed by X-ray and neutron diffraction. Although the 
distortion is suggested to be of Oj-type^^ for each Fe cube, 
the diffraction results themselves do not exclude a possibil- 
ity that the actural distortion is of a breathing type as found in 
PrRu4Pi2. We consider a situation where conduction electrons 
cause quantum fluctuations among degenerate AFQ configu- 
ration. If the macroscopic degeneracy is removed by this fluc- 
tuation mechanism, the cubic symmetry is kept on the av- 
erage, but a strong two-site AFQ correlation dominates the 
ground state. In order to pursue the idea qualitatively, we con- 
sider a toy model in which a conduction-electron spin is intro- 
duced at each site, and each spin has exchange coupling with 
a triplet. Namely the local interaction H^^l is given by eq.(7) 
for site A, and //gf^ for site B. Let us consider the system 

^r3+exc - Hr3 + ^exc + ^exc> (17) 

where we neglect the Fi singlet for simphcity. According to 
microscopic theoretical result,''' we take It < 0. In the Umit 
g3 = 0, we obtain the quartet ground state and the doublet ex- 
cited state for each site with separation 3|/t|/2. Thus the whole 
system is 16-fold degerate in the ground state. The AFQ in- 
teraction gj breaks the degeneracy. In the case of large g3/\It\, 
on the other hand, the lower manifold of 6 x 4 - 24 states 
has the following level sequence from the ground state: (i) 
quartet, (ii) octet, (iii) octet, (iv) quartet. The degeneracy of 
(i) comes from the spin degrees of freedom. Since there is 
no magnetic interaction between the pair, the Kramers de- 
generacy at each site cannot be broken. Figure 2 shows the 
energy levels as a function of gs/m taking the triplet as 
The asymptotes for \gi/It\ » 1 are characterized by the ferro- 
orbital state with the slope 784/3 as obtained earlier, and the 
antiferro-orbital state with the slope -196/ V3. The ratio of 
the slopes is -4/ V3 ~ -2.3. 

Now we include an intersite magnetic interaction of the 
form 
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Fig. 2. Energy levels of the triplet parr each of which is coupled with a spin 
in the conduction level. 

The ground-state quartet obtained above is split by this inter- 
action. A singlet ground state emerges in the case of /c < 0> 
and a triplet ground state otherwise. Figure 3 shows the lowest 
two levels for both signs of g^jlx. 

The wave function of the ground state for the system 
^^pair = Hr3+exc + Hc2 Can be derived exphcitly. We introduce 
notations: 

\s,;0} = 2-''H\'[}Mi}B + \i)A\l}B], (20) 

for a pair of CEF states and a pair of conduction spins, re- 
spectively. By cycUc rotation, we define similar states such as 
\yz)± and [s^; 0). The latter is given by 

\s,;0) = 2-"H\'[)A\'[)B + \i)A\i)B], (21) 

which has spin 1 with component Sx = 0. Then the ground 
state iff() of //pair with ^3 > 0, /t < and /c < is given by 

lAo = 3-1/2 0) + \y.)^\s^- 0) + \zxUsy; 0)] , (22) 

which does not depend on the interaction parameters because 
of its high symmetry. Although the pair has a complete AFQ 
correlation, there is no finite static quadrupole at each site. As 
a result, the cubic symmetry at each site is preserved on the 
average. 
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Fig. 3. The lowest two energy levels of the pair with inclusion of intersite 
coupling Jc of conduction spins. The numbers 1, 3 show the degeneracy of 
each level. 



Ha = J.s^^-sf\ 



(18) 



J. Phys. Soc. Jpn. 



Full Paper 



KuRAMOTO, Kiss, OxsuKiand Kusunose 5 



3.3 Staggered moments induced by magnetic field 

The dominance of quadrupolar moments in the ordered 
phase in PrFe4Pi2 is indicated by the resuhs of polarized neu- 
tron scattering^ where the induced AFM moments is about 2 
times larger for the field direction (001) than for (110). How- 
ever, there is no component of staggered magnetization per- 
pendicular to the magnetic field. If there is a static quadrupole 
order, one may naturally expect some anisotropy leading to 
the perpendicular component. The absence is not due to can- 
cellation by the mixture of quadrupole domains, since NMR 
has observed a sharp induced staggered field.^ 

If we assume the fluctuating state with no static AFQ or- 
der as described by the pair, there is no perpendicular com- 
ponent in the staggered magnetiztion. In our scenario of fluc- 
tuationg AFQ order, the static order parameter in PrFe4Pi2 
is a monopole, hexadecapole (rank 4 tensor) and/or a hexa- 
contatetrapole (rank 6 tensor). Both tensors can form a scalar 
in the poing group Th as given by O4, Og and Og in eq.(l). 
Similar multipoles have been postulated as order parameters 
in PrRu4Pi2. An important difl'erence, however, is the large 
Neel-like anomaly in the magnetic susceptibiUty in PrFe4Pi2, 
and its absence in PrRu4Pi2. We suspect the AFQ interaction 
is insignificant in PrRu4Pi2. 

3.4 Uniaxial pressure effect 

The quadrupole fluctuation is much affected by uniaxial 
stress, which serves as valuable probe to identify the nature 
of electronic states. According to ref.23, the susceptibility 
is much enhanced by uniaxial pressure along the magnetic 
field, while it is insensitive perpendicular to the field. The 
large anisotropy is common to both ordered and disordered 
phases. In order to understand the origin of such anisotropy, 
we have checked how the single- site triplet is affected by uni- 
axial stress. Provided the uniaxial stress along the z-axis fa- 
vors oblate charge distribution of 4/ states, the CEF triplet 
splits into the ground-state doublet and a singlet. Since the 
doublet |±) has a moment along the z-axis, the susceptibil- 
ity Xz is enhanced. On the other hand, Xx and Xy come from 
transition from the doublet and the singlet, which is reduced 
by the stress. Figure 4 shows an example of model calcula- 
tion for a Fi - r4 pseudo-quartet. Here the magnetoelastic 
couphng is taken to be gn - lOOK, and the elatic constant 
Cii - C12 = 35GPa.^^ It is seen that the decrease of the sus- 
cepibility for the perpendicular field is half of the increase 
along the field. This relation comes from the invariance of 
the average susceptibility in the present model. Although the 
tendency of anisotropy is consistent with experimental obser- 
vation, the quantitative discrepancy is serious. It seems nec- 
essary to include fluctuation effects even in the disordered 
phase. 

4. X-ray detection of hexadecapole order 

We discuss how to probe possible hexadecapole order tran- 
sition in PrFe4Pi2 and PrRu4Pi2. The most direct observation 
should be azimuthal angle scan in resonant X-ray scattering. 
In the quadrupole scattering process called E2, the staggered 
order of hexadecapoles contributes to the superlattice scat- 
tering. On the other hand, the hexacontatetrapole cannot be 
probed by the E2 scattering because of its higher rank (= 6). 
The combination O'^ + 50'^ constitutes a scalar in the cubic 
synnmetry, which has a staggered component in the hexade- 
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Fig. 4. Magnetic susceptibility of an isolated pseudo-quartet under uniaxial 
pressure cr = 3kbar. The parameters are taken as = 0.3K and A = 3K. 
The middle line shows the susceptibility without uniaxial stress. The unit 
of AT is (gjHsf/To with gj = 4/5 and To = 1 K. 
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Fig. 5. Predicted azimuthal angle dependence from the hexadecapole order 
in PrRu4Pi2 and PrFe4Pi2. Assumed are the (001) surface and the ct-tt' 
scattering. 

capole order. We have shown^ that the amplitude should de- 
pend on the azimuthal angle t// as cos 4i/' in the cr-cr' chan- 
nel, and as sin 4iJ/ in the a-n' channel. In the cr-n' channel, 
the polarization of the light is rotated by 90° from the inci- 
dent beam. These rapid oscillations of the amplitude are not 
shared by amplitude for other multipoles. In particular, there 
is no contribution from in the cr-n' channel. The result- 
ing intensity has the eight-fold pattern (oc 1 - cos 81^) with 
a minimum at t/' = 0. The cr-cr' channel will show a four- 
fold pattern by combination with O^. The intensity pattern is 
illustated in Fig. 5. In this way resonant X-ray scattering us- 
ing the E2 transition can identify the hexadecapole order. On 
the other hand, the non-resonant scattering also has contri- 
bution from hexadecapoles with angular momentum L - 4, 
which can be identified by the form factor analysis. Although 
the monopole contribution should dominate the intensity, one 
may hope to identify the contribution from L = 4 since the 
quadrupole component L = 2 is absent. Identification of the 
L = 4 component provides alternative detection of the hex- 
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Fig. 6. Comparison of the Schottky-type model for the specific heat with 
experimental results in ref.8. The unit of C is J/(mol K). 

adecapole order. 

5. Possible origin of heavy mass in PrOs4Sbi2 

The specific heat of PrOs4Sbi2 has a conspicuous anomaly 
of the Schottky type, which has been ascribed to the singlet- 
triplet CEF splitting. The superconducting transition occurs 
near the peak of this anomaly with a large jump (~ IJ/mol-K), 
which indicates a huge mass of the Cooper pair. If the su- 
perconductivity is destroyed by magnetic field, however, the 
heavy itinerant electrons seem to disappear at lower tem- 
peratures. Namely, the specific heat does not have the jT 
term with large y. Figure 6 compares the specific heat of the 
singlet-triplet CEF model in magnetic field with experimen- 
tal results.^ The field of // = 2T already masks the structure 
from the superconducting transition. It is seen that the dis- 
crepancy from the CEF model becomes larger for T > 1.5K, 
and the discrepancy above 4 K grows linearly with T, which 
gives the proportionality constant y ~ 0.9J/(mol-K^). In the 
low temperature end, on the other hand, the Schottky model 
gives reasonable account of the observed result. We there- 
fore conclude that a fraction of the entropy above ~1K is 
transferred to the temperature range higher than the maxi- 
mum of C(T). On the basis of observation above, we spec- 
ulate that a quadrupole Kondo efi'ect is active in this system 
where intra-triplet quadrupole excitations scatter conduction 
electrons. The mass enhancement from the intra-triplet exci- 
tation channel should naturally decrease if the thermal pop- 
ulation of the triplet becomes small. On the other hand, the 
singlet-triplet quadrupole channel should saturate below the 
temperature corresponding to the CEF sphtting. The lowest 
order efl'ect of this latter channel has been discussed by Fulde 
and Jensen.^* Since the yT term disappears at low temperature 
in Fig. 6, it seems that the intra-triplet channel is the dominant 
source of the heavy-fermion behavior. Microscopic study of 
^quad wiU shed further Ught on the vaUdity of this 

hypothesis. 

In contrast with the ordinary Kondo-type mechanism, the 
quadrupole scattering of conduction electrons with the 
symmetry requires a coherence of Pr atoms of the order of 
10 as shown by eq.(l 1). Hence it should be more sensitive to 
disorder such as La doping^^ than the ordinary Kondo effect. 
It should be a difficult but interesting problem how the heavy 



Cooper pairs survive or disappear in PrOs4Sbi2 at very low 
temperature where thermal population of the triplet becomes 
neghgible. 

6. Summary and Outlook 

In this paper we have proposed some ideas for understand- 
ing the intriguing properties of Pr skutterudites. The key to the 
understanding is the nature of the triplet CEF states and the 
singlet-triplet splitting which vary among PrFe4Pi2, PrRu4P|2 
and PrOs4Sbi2. The diff'erence of CEF states explains the 
presence of the strong Kondo effect in PrFe4Pi2, and the sharp 
CEF excitations observed in PrRu4Pi2 and PrOs4Sbi2. Re- 
cently Sm skutterudites have also been actively investigated. 
In particular, SmOs4Sbi2 has a large yT term in the specific 
heat, which is insensitive to magnetic field. We suspect that 
quadrupolar Kondo eff'ect may be the origin of the heavy 
mass. In contrast with the case of PrOs4Sbi2, the quadrupole 
active states are in the ground CEF level. These ideas can be 
corroborated only by microscopic calculation based on a real- 
isitic model. We plan to test the above ideas in the near future. 
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